Topology of the Fermi Surface Beyond the Quantum Critical Point 
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We examine the nature of phase transitions occurring in strongly correlated Fermi systems at the 
quantum critical point (QCP) associated with a divergent effective mass. Conventional scenarios for 
the QCP involving collective degrees of freedom are shown to have serious shortcomings. Working 
within the original Landau quasiparticle picture, we propose an alternative topological scenario for 
the QCP, in systems that obey standard Fermi liquid (FL) theory in advance of the QCP. Applying 
the technique of Poincare mapping, we analyze the sequence of iterative maps generated by the 
Landau equation for the single-particle spectrum at zero temperature. It is demonstrated that the 
Fermi surface is subject to rearrangement beyond the QCP. If the sequence of maps converges, 
a multi-connected Fermi surface is formed. If it fails to converge, the Fermi surface swells into 
a volume that provides a measure of entropy associated with formation of an exceptional state 
of the system characterized by partial occupation of single-particle states and dispersion of their 
spectrum proportional to temperature. Based on this dual scenario, the thermodynamics of Fermi 
systems beyond the QCP exhibits striking departures from the predictions of standard FL theory. 
Mechanisms for the release of the entropy excess of the exceptional state are discussed. 

PACS numbers: 71.10.Hf, 71.27. +a, 71. 10. Ay 



I. INTRODUCTION 

The Landau quasiparticle patternA of low-temperature 
phenomena in Fermi systems is universally recognized 
as a cornerstone of condensed-matter theory. Although 
Landau's Fermi liquid (FL) theory was originally formu- 
lated for liquid 3 He, it is a quirk of fate that discrep- 
ancies between theoretical predictions of the theory and 
experimental data on the low-temperature behavior of 
key thermodynamic properties first came to light in this 
system. 

In the three-dimensional (3D) liquid, deviations of the 
spin susceptibility x(T) and specific-heat ratio C{T)/T 
from the predicted constant behavior, rather small at 
very low temperatures T, increase with T. In 2D liq- 
uid 3 He, however, departures from the predictions of FL 
theory become more pronounced as the temperature is 
lowered once the density p reaches a critical region where 
the FL effective mass M* (p) is greatly enhanced* 2 -^ the 
values of M* having been extracted from the specific heat 
data via the FL formula C FL (T) = Tp F M*/3. This con- 
trary behavior rules out damping of single-particle exci- 
tations as the cause of the observed failure of FL the- 
ory in 2D liquid 3 He, since damping effects decrease and 
vanish at T — » 0. The experimental data on the spin 
susceptibility x{T) of 3 He films present a thorny chal- 
lenge of interpretation. In agreement with the Curie law 
but against the FL predictions, the product x(T)T fails 
to show vanishing behavior as very low values of T are 
reached in the measurements. In the critical density re- 
gion, the value of this product gradually increases with 



increasing p, as if a fraction of spins of 3 He atoms be- 
comes localized and coexists with the liquid part of the 
system. Analogous behavior has been observed for 2D 
electron gas£ Non-Fermi-liquid (NFL) behavior also ex- 
hibits itself in properties of strongly correlated electron 
systems of heavy-fermion metals ^^ ' 10 i 11 ' 12 i 13 ' 14 ' 15 

Various theories, generally invoking critical fluctua- 
tions associated with second-order phase transitions ; 1 ^ 17 
have been proposed to explain NFL behavior in strongly 
correlated Fermi systems beyond the quantum critical 
point (QCP). At a QCP, which by definition occurs at 
zero temperature, the Landau state becomes unstable 
due to a divergence of the effective mass M * (cf. Ref. [Lsl 
and works cited therein). 

On the weakly correlated side of the QCP, the proper- 
ties of systems of interest are described within the stan- 
dard FL theory, in which the Fermi liquid is treated 
as a gas of interacting quasiparticles with a spectrum 
of single-particle excitations given by e(p) — pf(p — 
Pf) /M* . Since this spectrum loses its meaning when M* 
diverges, the standard FL theory fails at the QCP and 
beyond. Bearing in mind the fundamental role played 
by the Landau approach in modern condensed matter 
physics, there is ample incentive to investigate the situa- 
tion beyond the QCP within the original Landau quasi- 
particle formalism ji free of shortcomings of the standard 
quasiparticle picture. We will demonstrate that if the 
strength of the effective interaction between quasiparti- 
cles reaches a critical level, a topological phase transition 
occurs, and the properties of the system change dramat- 
ically. 

This paper is organized as follows. In Sec. II we call 
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attention to certain basic flaws of the conventional collec- 
tive scenario for the QCP that envisions a "fatal" break- 
down of the quasiparticle picture when the quasiparticle 
weight z vanishes at a second-order phase transition. We 
then describe an alternative topological scenario for the 
QCP, in which the quasiparticle group velocity at the 
Fermi surface changes sign while the z -factor stays fi- 
nite. In Sec. Ill we apply the concept and technique of 
Poincare mapping, widely used in the analysis of non- 
linear phenomena, to explore the structure of approxi- 
mants (iterates) generated by iteration of the nonlinear 
integral equation^ for the single-particle spectrum e(j>) 
and momentum distribution n(p) at zero temperature. 
Beyond the QCP, the iteration procedure either con- 
verges, in which case the Fermi surfaces becomes multi- 
connected; or it does not. In the latter case, which applies 
for systems with long-range effective forces, the patterns 
of successive iterates for e(j>) and n{p) acquire chaotic 
features. A special procedure is introduced for averag- 
ing over sequences of iterates, in such a way as to sup- 
press these chaotic features. In Sec. IV, we demonstrate 
that the resulting averaged single-particle energies and 
occupation numbers coincide with those belonging to a 
state with a fermion condensate-an exceptional type of 
ground state possessing a distinctive topological struc- 
ture. Sec. V presents results from numerical calculations 
based on the Landau equation for the spectrum e(p) at 
finite temperatures beyond the points of the topological 
phase transitions. The relevance of the topological sce- 
nario to the real experimental situation in Fermi liquids 
beyond the QCP is discussed in Sec. VI. Extraordinary 
aspects of the NFL thermodynamics of these systems — 
Curie- Weiss behavior of the spin susceptibility and the 
role of the entropy excess associated with the fermion 
condensate — are addressed in Sec. VII. In Sec. VII, we 
analyze different possibilities for release of this excess 
entropy, which may be responsible for the diversity of 
quantum phase transitions appearing in the phase dia- 
grams of strongly correlated Fermi systems beyond the 
QCP. The paper is concluded in Sec. IX with a summary 
of key developments and findings, and with remarks on 
the true scope of Landau theory. 



II. TWO DIFFERENT SCENARIOS FOR THE 
QUANTUM CRITICAL POINT 

A dominant activity in condensed-matter physics dur- 
ing the last decade has been the investigation of quantum 
phase transitions, occurring at extremely low tempera- 
tures in strongly correlated Fermi systems. Imposition 
of pressure or magnetic fields allows one to push the 
transition temperature toward zero, producing a QCP 
associated with divergence of the density of states, or 
equivalently, the effective mass M*{p) at a critical den- 
sity Poo- Li homogeneous nonsuperffuid Fermi systems, 
the ratio of M* to the bare mass M is determined by the 



textbook formula 
M 



(1) 



where = p 2 /2M — /i is the bare quasiparticle energy 
measured relative to the chemical potential, S repre- 
sents the mass operator, and the quasiparticle weight 
z in the single-particle state is given by z — [1 — 
(9E(p, e)/de) Q ]~ 1 . Here and henceforth, the subscript 
indicates that the quantity in question is evaluated at 
the Fermi surface. 

Since the Feynman-Dyson era, it has been a truism 
that the effects of the two factors in Eq. |T]) , associated re- 
spectively with (i) the energy dependence of the mass op- 
erator (or self-energy) X(p, e) and (ii) its momentum de- 
pendence, cannot be separated from each other based only 
on measurements of thermodynamic and transport prop- 
erties. Importantly, at the quantum critical point these 
factors express themselves differently in different scenar- 
ios. In a conventional collective scenario for the QCP, the 
energy dependence of S plays a decisive role. "Quasipar- 
ticles get heavy and die,"-iS since critical fluctuations de- 
stroy the quasiparticle picture, causing the quasiparticle 
weight z to vanish at the transition point i 19 ' 20 i 21 

By contrast, in a topological scenario 2 ^, for the QCP, 
which is associated with a change of sign of the quasi- 
particle group velocity vf = (de(p)/dp)o appearing in 
Eq. ((T|), the momentum dependence of the mass opera- 
tor evidently assumes a key role, when we note that Vf 
is proportional to the sum 1 + (9S(p, e = 0)/c?e°) Q . In 
this scenario, nothing catastrophic happens beyond the 
critical point where vf reverses sign: the original quasi- 
particle picture, advanced by Landau in the first article 
devoted to theory of Fermi liquids^ holds on both sides 
of the QCP, since the z-factor remains finite. 



A. An original Landau quasiparticle pattern 

We recall that the heart of the Landau quasiparticle 
picture is the postulate that there exists a one-to-one cor- 
respondence between the totality of real, decaying single- 
particle excitations of the actual Fermi liquid and a sys- 
tem of immortal interacting quasiparticles. Two features 
specify the latter system. First, the number of quasi- 
particles is equal to the given number of particles (the 
so-called Landau-Luttinger theorem). This condition is 
expressed as 



n(p) dv = p , 



(2) 



where n(p) is the quasiparticle momentum distribution, 
the factor 2 comes from summation over the spin projec- 
tions, and dv is a volume element in momentum space. 
Second, the entropy S, given by the combinatorial ex- 
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pression 



S(T) 



[n(p) lnn(p) + (1— n(p)) hi(l— n(p))] dv , 

(3) 

based on the quasiparticle picture, coincides with the en- 
tropy of the actual system. Treating the ground-state 
energy Eq as a functional of n(p), Landau derived the 
formula 



i(p,T) = 



1 



z(p)/T 



(4) 



which has an obvious (but misleading) resemblance to 
the Fermi-Dirac formula for the momentum distribution 
of an ideal Fermi gas. In contrast to the ideal-gas case, 
the quasiparticle energy e(p) = SEo/Sn(p) — p itself must 
be treated as a functional of n(p). 

Another fundamental relation, which follows from the 
Galilean invariance of the Hamiltonian of the problem, 
connects the group velocity de(p)/dp of the quasiparti- 
cles to their momentum distribution through ^ 23 ' 24 



de(p) 
dp 



P_ 

M 



1 [ft \ dn (Pl) A 



(5) 



The interaction function /(p,Pi) appearing in this rela- 
tion is the product of z 2 and the scalar part of the scatter- 
ing amplitude T u . In turn, T u is the w-limit of the scat- 
tering amplitude T of two particles whose energies and 
incoming momenta pi,P2 lie on the Fermi surface, with 
scattering angle cos6> = PiPilvp and the 4-momentum 
transfer (q, uj) approaching zero such that q/ui — > 0. 

Landau (see formula (4) in Ref. 1) supposed that so- 
lutions of Eq. ([5]) always arrange themselves in such a 
way that at T = 0, the quasiparticle group velocity vf 
maintains a positive value, implying that the quasipar- 
ticle momentum distribution takes the Fermi-step form 
n{p,T = 0) = np(p) = 0{p — pp). If this supposition 
holds, implementation of the original quasiparticle pic- 
ture is greatly facilitated: properties of any Fermi liquid 
coincide with those of a gas of interacting quasiparticlcs. 

The failure of this assumption in strongly correlated 
Fermi systems, established first in microscopic calcu- 
lations of Refs. I25ll26l . can be seen from the analysis 
of Eq. ^ itself. Consider for example homogeneous 
fermionic matter in 3D. Upon setting T = and p = pp 
in Eq. (O and denoting the first harmonic of the inter- 
action function by /1 we find 



vf 



P_F 

M 



(6) 



having introduced the dimensionless parameter F® — 
PfM/i/tt 2 . Eq. ([6]) is easily rewritten in the FL 
formic 



M 
~M* 



lp F M 



= 1 - - 

3 7T 



1 



fl = 1 - o F l 



(7) 



Hence the inequality vp > is violated at the critical 
density where -F^Poo) = 3. 



It is instructive to rewrite Eq. ([7|) in terms of the de- 
limit of the dimensionless scattering amplitude vT k — 
A + Series, where v — z 2 pfM* /tt 2 is the quasiparti- 
cle density of states. Noting the connection 1 ^ 23 ' 24 A\ = 
Fi/(l + Fi/3) oiA x to Fi =p F M*fi/n 2 , simple algebra 
based on Eq. ([7]) then yields 



M 1 1 A 

= 1 Ai . 

M* 3 

Clearly then, one must have 

^i(Poo) = F 1 (p oo ) = 3 



(8) 



(9) 



at the density where the effective mass diverges. 

We would like to emphasize that the critical density 
has no relation to the density associated with violation of 
the Pomeranchuk condition for stability against dipolar 
deformation of the Fermi surface. In the latter situation, 
the critical value M* c of the effective mass M* is found 
from the condition 1,23,24 



VfM{ c 
; 3tt 2 







(10) 



which, according to Eq. (O, may be recast as (1 — 
F^/3)" 1 = M*JM = 0. Quite evidently, this situation 
is not relevant to the QCP. 



B. Critique of the conventional QCP scenario 

In this subsection, we show that the collective scenario 
for the QCP encounters difficulties when the wave vector 
k c specifying the spectrum of the critical fluctuations has 
a nonzero value. This result casts doubt on the standard 
collective QCP scenario, since the assumption of a finite 
critical wave vector number k c is aligned prevalent ideas. 



1. Homogeneous matter 

We first address the homogeneous case, and, for the 
sake of specificity, we focus on density fluctuations 2 ^ as 
the seed for the transition. The requirement of antisym- 
metry of the amplitude T with respect to interchange of 
the momenta and spins of the colliding particles leads to 
the relation 2 ^ 

A( Pl ,p 2) k,a; = 0;p^p c ) = -D(k) + ^D(p 1 -p 2 + k) , 



where 

D(k — > k c , lo = 0; p — * p c ) 



t- 2 {p) + {k~k c ) 2 



(11) 
, (12) 



with g > and the correlation length £(ft) divergent at 

P = Pc- 
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To explicate difficulties encountered by the standard 
scenario for the QCP in homogeneous matter, let us cal- 
culate harmonics Ak(p) of the amplitude A(pp,pp, cos 9) 
from Eqs. (fTTj) and (jT2J) . In particular, we obtain 



A (p^p c ) 



7T k c £,(p) 



Ai(p^Pc) 



3tt k c £(p) 

g— — cos 6> 

2 P F 



(13) 

The sign of Ai(p — > p c ), which coincides with the sign 
of cos^o = 1 — k~/2p F , turns out to be negative at 
k c > PfV%- According to Eq. (jSJ, this means that at 
the second-order phase transition, the ratio M*(p c )/M 
must be less than unity. We are then forced to conclude 
that the densities p c and poo cannot coincide. Although 
the z-factor does vanish at the density p c due to the diver- 
gence of the derivative (dT,(p,e)/de) , the effective mass 
M* remains finite, since the derivative (d£(p, e)/9e°) 
diverges at the QCP as well^ 

Moreover, it is readily demonstrated that at positive 
values of cos6>o below unity, vanishing of the z-factor is 
incompatible with divergence of M*. Indeed, as seen from 
Eq. (fT3"|) , the harmonics Ao(p c ) and Ai(p c ) are related to 
each other by A (p c ) — Ai(p c ) / (3 cos 9o) ■ IfM*(p c )were 
infinite, then according to Eq. |8[). A\{p c ) would equal 
3, while Aq(p c ) = l/cos^o > 1. On the other hand, 
the basic FL connection Aq — Fq/(1 + Fq) implies that 
Aq < 1, provided the Landau state is stable. Thus in 
the conventional scenario, the QCP cannot be reached 
without violating the stability conditions for the Landau 
state. The same is true for critical spin fluctuations with 
nonvanishing critical wave number. One must conclude 
that the system, originally obeying FL theory, undergoes 
a first-order phase transition upon approaching the QCP, 
as in the case of 3D liquid 3 He. 

Based on these considerations we infer that for homo- 
geneous matter, the list of possible second-order phase 
transitions compatible with the divergence of the effec- 
tive mass M* includes only long- wave transitions. These 
are associated with some L— deformation of the Fermi 
surface such that one of the two Pomeranchuk stability 
conditions^ e.g., 1 + F L /(2L + 1) > 0, is violated. 



2. Anisotropic systems 

We turn now to systems, typified by heavy-fermion 
metals, whose Fermi surface is anisotropic. In the en- 
ergy region |e| ~ T most relevant to the thermodynamic 
and transport properties of all Fermi systems, includ- 
ing anisotropic examples, the quasiparticle group velocity 
v = de/ dp consists only of the component v n normal to 
the Fermi surface where e = 0; hence it depends only on 
the corresponding momentum component p n . In homo- 
geneous matter, the value of the velocity is the same at 
any point of the Fermi surface, whereas in the anisotropic 
case, it depends on the observation point. Even so, the 
strategy used below to investigate topological phase tran- 
sitions in anisotropic systems is much like that applied 



above in the analysis of such transitions in homogeneous 
matter. 

To deal with electronic systems of solids, Eq. (fT|) is 
replaced by 



«n(p) 



de(p) 

dp n 



de°(p) 

dp n 



1 



/3£(p,e) 



V deO(p) , „ 

(14) 

where e°(p) is the "bare" electron spectrum, evaluated 
so as to account for the external field due to the lat- 
tice, while the electron mass operator S includes all the 
electron-electron interaction effects. 

Violation of Galilean invariance in solids results in a 
change of Eq. ([5]) stemming from the FL relation be- 
tween the k- and w-limits of the vertex T(p), because 
the Pitaevskii identity 2 ^ zT w (p) = p used in deriving 
this equation is no longer valid. Accordingly, we arrive 
at the generalized version 



dp 



<n Jo J Opi n 



where dS is the surface element in momentum space. In 
obtaining this relation, we have neglected a possible p 
dependence of the interaction function /; this simplifi- 
cation allows a straightforward integration over coordi- 
nates. Additionally, we have employed the identity^ 



C^(p))o 



( dG~\p) 
V dp 



de{p) 
dp 



(16) 



implied by the gauge invariance of the problem, where G 
is the single-particle Green function. 

Assuming the "bare" spectrum e°(p) to be stable, po- 
tential points of divergence of the density of states are 
associated — as in the preceding discussion of alternative 
scenarios for the QCP in homogeneous matter — with (i) a 
vanishing of the renormalization factor z and consequent 
breakdown of the quasiparticle picture or (ii) a change of 
sign of the group velocity v n , leading to a change of the 
connectivity of the Fermi surface. 

We restrict the analysis to the much-studied case 
of a quasi-two-dimensional heavy-fermion metal with a 
quadratic lattice. We also assume filling close to 1/2. 
Conventional belief would hold that in this case, antifer- 
romagnetic fluctuations play the key role in the structure 
of the corresponding QCP . 18 i 21 The relevant propagator 
takes the form 



D(k) 



1 



r 2 + (k-Q)^ 



(17) 



with the critical wave vector Q = (n / a,~K / d), where a 
is the lattice constant. However, this viewpoint is insup- 
portable. To some extent, its weakness is already evident 
from the above treatment of homogeneous matter; tak- 
ing account of fluctuations with k c > pfV% results in a 
suppression of the effective mass M*. Attending now to 
the anisotropic system, based on the findings of Ref. [13, 




FIG. 1: Location of roots of the equation e(p, T = 0) = mea- 
sured from the chemical potential p, beyond the bifurcation 
point 



it may be shown that in the postulated situation z — > 0, 
the first term on the right-hand side of Eq. (| 15(1 reduces 
to de°(p)/dp n to yield 



de(p) 

dp n 



de°(p) 

dp n 



gn(gi)/gpiw 

e- 2 + (P-Pi-Q) 



2 dSidpi n , 

where g > 0. In the case of ferromagnetic or long- wave 
antiferromagnetic critical fluctuations, the overwhelming 
contributions to the integral term of Eq. (fT5|) come from 
a region pi ~ p where the derivatives dn{p)/dp n and 
de°(p)/dp n have opposite signs. Consequently, account- 
ing for these critical fluctuations leads to a suppression of 
the group velocity v n , promoting emergence of the QCP. 

What happens if instead one pursues the case of short- 
wave antiferromagnetic fluctuations, invoked in conven- 
tional scenarios of the QCP? In this case, the so-called 
hot spots — points of the Fermi surface connected by 
the antiferromagnetic vector Q — contribute predomi- 
nantly to the integral term of Eq. (fT8|) . The quantities 
de a (p)/dp n and dn(p— Q)/dp n now have the same sign 
and thus cooperate to enhance rather than suppress the 
value of the group velocity v n (p) . We see then that an- 
tiferromagnetic fluctuations actually have the effect of 
hampering the emergence of the QCP in anisotropic sys- 
tems, in correspondence with the result obtained for ho- 
mogeneous matter. Accordingly, despite common belief, 
in heavy-fermion metals the QCP cannot be attributed 
to antiferromagnetic fluctuations with the critical wave 
vector Q = (ir/a,Tr/a). 



C. Salient features of the topological scenario for 
the QCP 

In ordinary ("canonical") Fermi liquids, there exists 
a one-to-one correspondence between the single-particle 
spectrum e(p) and the momentum p, at least in the vicin- 
ity of the Fermi surface. This is equivalent to asserting 
that the equation 
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FIG. 2: Single-particle spectrum e(p, T = 0) of a homoge- 
neous 2D electron gas, measured from the chemical potential 
p, in units of £%• = p%/2M, evaluated microscopically^ at 
T = for different values of the dimensionless parameter 
r s = V2Me 2 /p F - 



has a single solution p = pp, since the group velocity vf 
is positive. Within FL theory, Eq. (|19[) reduces to the 
relation 



e(p,T = 0;n F ) 







(20) 



E(p, e = 0) = p 



(19) 



where the energy e(p), measured from the chemical po- 
tential (i, is evaluated with the Landau quasiparticle mo- 
mentum distribution uf(p) — &(p — Pf)- However, it is 
a key ingredient of the topological scenario for the QCP 
that at a critical value of some input parameter, specifi- 
cally the density or a coupling constant 1700 , the group 
velocity vp vanishes. (This feature in demonstrated in 
several numerical examples presented in Sec. V.) 

Beyond the critical point, e.g. at g > g^, Eq. (j2"0|) 
acquires at least two new roots (see Fig. [1]), triggering 
a topological phase transition. 28 Significantly, terms pro- 
portional to elne, which are present in the mass operator 
E of marginal Fermi liquids, do not enter Eqs. (|19|) and 

The bifurcation point p^ of these equations can emerge 
anywhere in momentum space. Here we examine the 
case in which p^ coincides with the Fermi momentum 
Pf, so that at the critical density Poo, the group veloc- 
ity vp(poo) changes its sign. This distinctive feature of 
the topological scenario for the QCP is in agreement with 
results of microscopic calculations of the single-particle 
spectrum e(p, T = 0) of 2D electron gas^ as illustrated 
in Fig. [21 It is seen that the sign of vf remains positive 
until the dimensionless parameter r s attains a critical 
value Too ~ 7.0, where vf vanishes. At larger r s , the 
group velocity vf = (de(p,n,F)/dp) Q , evaluated with the 
quasiparticle momentum distribution np{p), is negative. 
Having lost it stability, the Landau state is replaced by a 
new state through the intervention of a topological phase 
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transition. 

A qualitatively similar situation is to be expected in 
anisotropic Fermi systems. In support of this statement, 
we may refer to the first example of the QCP, namely 
the Lifshitz saddle point, 3 - which is characterized by a 
density of states diverging at T — > as N(T, p^) ~ 
ln(l/T). This example has been elaborated specifically 
for anisotropic 2D electron systems in solids. Close to 
the saddle point, the QCP single-particle spectrum has 
the form e(p, Pco ) = {l/2)[p 2 JM xx - p 2 y /M vy ), with the 
components p x and p y specifying the distance to the sad- 
dle point in momentum space. In the Lifshitz model, 
the quantities M xx and M yy have comparable finite val- 
ues. If one of these parameters, say M xx , becomes much 
larger than M yy , we arrive at the so-called extended sad- 
dle point, which has been considered in connection with 
high-T c superconductivity in Ref. l3ll . The corresponding 
phase transition, associated with the divergence of M xx , 
is also topological. 



III. POINCARE MAPPING FOR STRONGLY 
CORRELATED FERMI SYSTEMS 



In this section we restrict considerations to homoge- 
neous matter at zero temperature. Microscopic calcula- 
tions of the single-particle spectrum e(p, T = 0) are as 
yet available only for several simple types of bare inter- 
actions between particles. Moreover, there is no micro- 
scopic theory beyond the point where the Landau state 
loses its stability. On the other hand, given a phenomeno- 
logical interaction function /, Eq. ([5]) holds on both the 
sides of the topological QCP, since the z-factor retains 
a nonzero value. Numerical iteration is a standard ap- 
proach to solution of an equation such as Eq. ([5]). We 
shall see that the mathematical counterpart of iteration, 
Poincare mapping, which has been widely exploited in 
theory of turbulence^ 2 , is also instrumental in elucidat- 
ing the striking features inherent in solutions of Eq. ([5]) 
beyond the QCP. 

The discrete iterative map corresponding to Eq. ([5]) 
reads 



de^ +1 \p) 
dp 



P_ 

M 



1 

3^ 



hip, Pi) 



dpi 



p\dpi 



(21) 

the iterate fj,^ +1 ^ for the chemical potential being deter- 
mined from the normalization condition ([2|). The index 
j = 0, 1, 2, . . . counts the iterations (zeroth, first, second, 
. . .). The iterate n^ +1 ^(p) of the momentum distribution 
n{p) is generated by inserting the corresponding spectral 
iterate e^ +1 \p) into Eq. (U|), which at T = reduces to 
a Heaviside function n(p) = 9(—e(p)). 



A. Poincare mapping for systems with a 
nonsingular interaction 

The key quantity of the Poincare analysis in the Fermi 
liquid problem is the group velocity 



de^jp) 
dp 



JL 

M 



3tt 2 



fi{p,PF\p) 



(22) 



evaluated by inserting the standard FL momentum dis- 
tribution riFip) — 9{p — pf) into the right-hand side of 
Eq. (|2~lj) as the zero iterate (starting approximant) for the 
quasiparticle momentum distribution nip). In canonical 
Fermi liquids, for which the sign of vp is positive and the 
function e^ip) has the single zero at p = pp, the first 
iterate (p) and all higher iterates for the distribution 
nip) coincide with the np(p) 7 this being a fixed point of 
the transformation. However, beyond the QCP the sign 
of the group velocity vp evaluated from Eq. ((22|) becomes 
negative, and the first iterate e^\p) = e(p;np) for the 

spectrum already has three zeroes p± < jxp < p!j , 
implying three kinks in the momentum distribution. 

At the second step, the evolution of the iteration pro- 
cess follows different patterns, depending on the presence 
or absence of a long-range component in the effective in- 
teraction between quasiparticles (long range in coordi- 
nate space). We first consider the simpler case in which 
the effective interaction /(fc), local in coordinate space, 
has no singularities. In this case the interaction function 
fi(PiPF) can then expanded in a Taylor series in the 
variable p — Pf, and Eq. §5§ may be recast as a set of al- 
gebraic equations. Then, as the straightforward analysis 
demonstrates and the numerical calculations described in 
Sect. V confirm, the shape 



e(p,T= 0) ~ (p - pt)(p - p 2 )ip - Pa) 



(23) 



of the spectrum remains the same independently of the 
number of iterations, with three roots p\ < P2 < P3 spec- 
ifying the location of the three sheets of the Fermi surface 
and lying close to the QCP Fermi momentum p^. Cor- 
respondingly, the spectrum e(p, T — 0) changes smoothly 
in the momentum regimes removed from the kinks, but 
oscillates in the interval [pi,j>3]. The amplitude of the 
oscillation, i.e., the maximum value 

T m = max \e(p, T = 0)\, p x < p < p 3 , (24) 

of the departure of \e(p,T — 0)| from 0, defines a new 
energy scale of the problem: if temperature T attains 
values comparable to T mi the kink structure associated 
with the multi-connected Fermi surface is destroyed. 

To evaluate T m , we note that according to Eq. (|2"3"|) 
the group velocity de(p)/dp is a parabolic function of p, 
conveniently written as 



dejp) 
dp 



-X + Y(p-p 



2 

m ) i 



(25) 
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where p m determines the location of the minimum of the 
group velocity. Comparison of Eqs. (j2"3"j) and (j2"5|) leads 
to the following set of equations, 



n + r 2 + r 3 = , 



nr 2 + nr 3 + r 2 r 3 



3A 



(pm + nf - (p m + r 2 f + (p m + r 3 ) 3 = 37r 2 p , (26) 

for the three (shifted) roots rk = Pk—p-m, with k = 1, 2, 3. 
The last of Eqs. (|26p is obtained from the normalization 
condition ((2|). 

In the vicinity of the QCP where X = X^. = 0, the 
parameter p m and the positive quantities X(p) and Y(p) 
can be shown to change linearly with p — p^ . Simple but 
lengthy algebra then yields 



r 3 



-r\ ~ I/O - Poo\ 1/2 , r 2 ~\p-poo\ 



(27) 



These results allow us to express relevant parameters in 
terms of the difference p — p^ , namely (i) Ap — p 3 — pi , 
which characterizes the range of the flattening of the 
spectrum e(p, T = 0) beyond the QCP, (ii) the tempera- 
ture T m associated with the crossover from standard FL 
behavior to NFL behavior, and (iii) the zero-temperature 
density of states iV(0) ~ Sfc[^ e (P = P^/dp] -1 , which re- 
places the ratio M* /M in the standard FL formulas for 
the specific heat, spin susceptibility, etc. We find 



Ap 



\P~Po 



11/2 rp 



\p- Poo \ 3 / 2 , N(0) 



\P-Poc\ 1 ■ 

(28) 

To our knowledge, the multi-connected Fermi surface, 
shown here to arise in the homogeneous case, was first 
considered in 1950 by Frohlich. 33 Within the Hartree- 
Fock (HF) approach, model variational calculations for 
the ground-state energies of homogeneous systems lead- 
ing to a multi-connected Fermi surface were performed 
more than 20 years ago^ while the corresponding phase 
transition was first discussed in terms of HF single- 
particle spectra in Ref. [3^ and later in Ref. [36i. The 
original calculations of properties of this transition on 
the basis of Eq. ([5]) were performed in Refs. [37ll38l . These 
calculations show that as the coupling constant moves 
away from the QCP value g^, the number of sheets of 
the Fermi surface, which coincides with the number of 
the roots of Eq. (|20|) , grows very rapidly, with the dis- 
tance between the sheets shrinking apace. Importantly, 
in all these solutions, the relation 

n 2 {p) = nip) , (29) 

inherent in the standard FL picture, is still obeyed. 



B. 2-cycles in Poincare mapping for systems with a 
singular interaction function / 

A remarkable feature of equations relevant to the tur- 
bulence problem is the doubling of periods of motion, 



P'P F 

0.9 1.0 1.1 
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1.0 



CL 
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FIG. 3: Iterative maps for the quark-gluon plasma problem 
(|3ip with the bare spectrum e° = cp and the dimensionless pa- 
rameter a = g/c = 0.3. Left panels: spectral iterates e' J '(p) 
with j — 0,1,2, in units of cpF- Right panels: momentum- 
distribution iterates n"'(p). The 2-cycle reveals itself in the 
coincidence between the first and third right panels. 



giving rise to dynamical chaos . 32 ' 39 If now the nonlinear 
system corresponding to Eq. (f2Tj) is considered within 
this context by associating an iteration step with a step 
in time, one's first instinct is to assert that such a phe- 
nomenon cannot occur when Poincare mapping is imple- 
mented, since chaos in the classical sense cannot play a 
role in the ground states of Fermi liquids at T = 0. And 
indeed, this general assertion seems to be validated by the 
results of the preceding subsection. However, the Tay- 
lor expansion of the interaction function ,Ji{p,pf) that 
we employed in the above analysis fails in the case of 
long-range effective interactions. The Fourier transform 
f(k — > 0) becomes singular, and the previous analysis is 
inapplicable. 

As an interesting example, we note that such a singu- 
larity exists in the bare interaction 



Pi ■ P2 " (Pi ■ k)(p 2 ■ k)/fc 2 
fc 2 



r°(p 1 ,p 2 ,k,w = o) = -5 

(30) 

between quarks in dense quark-gluon plasma, wherein 
g > 0. In this example, Eq. ([5]) has the fori 



,35 



de(p) _ de a p 



g / In 



2pp dn(p 



-dpi , (31) 



dp dp J p — pi dpi 

where e° ~ cp is the bare single-particle spectrum of light 
quarks. 
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The first iterate for the spectrum, evaluated from 
Eq. (f3"Tj) with the distribution np(p), has an infinite neg- 
ative derivative de^ /dp at the Fermi surface. It then 
follows that independently of the value of the dimension- 
less coupling constant a — g{de^/ dp)^ , Eq. (120]) . with 
x = p/pF — 1) takes the form 



x 1 



a In - ) =0 

x . 



(32) 



and has three different roots — xq, 0, and xq, with xq = 
2 exp -1 /". The corresponding first iterate of the momen- 
tum distribution is n^(p) = 9(x + x ) — 8{x) +9(x — x ). 
The next iteration step yields n^(p) = rip(p), so the 
standard FL structure of the momentum distribution is 
recovered. The nonlinear system enters a 2-cycle that 
is repeated indefinitely; one is dealing with a 2-cycle 
Poincare mapping. 

The first two iterations of the mapping process are 
illustrated in Fig. [3] The top-left panel of this figure 
shows the bare spectrum e^(p) = e°(p). The first iter- 
ate of the spectrum, e^(p), appearing in the middle- left 
panel, is evaluated by folding the kernel \n(2pp/ (p — p\j) 
with the Fermi-step Hf(p), shown in the top-right panel. 
The spectrum e^(p) possesses three zeroes: p\ < pf, 
Pi = Pf, and p% > pp, implying that the first iterate 
n^\p), drawn in middle-right panel, describes a Fermi 
surface with three sheets. This distribution differs from 
the ordinary Fermi step only in the momentum interval 
— xq < x < xq. The next iterate, e^(p) (bottom-left 
panel), again has a single zero pf, and the correspond- 
ing momentum distribution (p) (bottom-right panel) 
coincides identically with np{p). 

A 2-cycle Poincare mapping also arises in treating the 
Nozieres model; 41 - for which the interaction function / 
has the limited singular form /(q) = gS(q) with g > 0. 
In this model, the iterative maps (illustrated in Fig. 
are generated from the equation 

(p) + = P 2 /2M + gn® (p) , (33) 



along with the normalization condition ^ for 
n 0'+i)(p) = 0(_ e O'+i)( p )). Here, the odd iterates 
n ( 2 3+ 1 ) (p) of the momentum distribution deviate from 
the np(p) in the interval —g/4e° F < x < g/4s F , but in 
even iterations, the Fermi step reappears intact. 

Numerical analysis demonstrates that similar 2-cycles 
arise when Poincare mapping based on Eq. (|2ip is imple- 
mented for other systems having an interaction function 
that is singular at k — > (see Sec. V). In all these cases, 
the emergence of a 2-cycle turns out to be an unam- 
biguous signal of the instability of the standard Landau 
state. Moreover, the domain of momentum involved in 
the cyclic behavior is almost identical with the domain 
within which an improved iteration algorithm fails to con- 
verge (see Sec. III.D). As will be seen, this concurrence is 
significant in that the associated volume of momentum 
space provides a measure of the entropy stored in the 
exceptional ground state that replaces the Landau state. 
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FIG. 4: Same as in Fig. [3] for the discretized Nozieres model 



33]) with g = 0.2e u F . 



C. 2-cycles in Poincare mapping for finite systems 



In finite Fermi systems — nuclei, atoms, atomic clus- 
ters, quantum dots, molecules, . . . — integration over mo- 
menta in Eq. (|2ip is replaced by summation over single- 
particle quantum numbers. As a concrete illustration, 
let us explore a schematic model^ of a spherical nucleus, 
in which the single-(quasi)particle energies are indepen- 
dent of the magnetic quantum number m. We consider 
two neutron levels available for filling in an open shell, 
denoted and +. 

In the presence of a quasiparticle interaction /, the en- 
ergies of the two levels arc influenced when neutrons are 
added to the system. With this in mind, the energies 
of levels and + when N neutrons are added to level 
and none added to level + are denoted €0(^,0) = cq(N) 
and €+{N, 0) = e + (N), respectively. The initial dis- 
tance between the two levels, with no added neutrons, 
is D{0,0) = e+(0,0) - e (0,0) = D > 0, having des- 
ignated as the lower of the two levels. The behav- 
ior of the level distance D{N,0) = e + (N,0) - e (AT,0) 
as neutrons are added is crucial to the behavior of the 
system. To facilitate the argument, we adopt a highly 
simplified interaction between neutrons, retaining in / 
only a principa l 43 ! 44 5(r)-like component. Also for sim- 
plicity, the relevant matrix elements for determining the 
energy shifts are reduced to two, namely /oo = /++ and 
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fo+ = f+o, which (with k = 0, +) are calculated as 

fk*= J Rl{r)f{p{r))Rl{r) r ^- (34) 

in terms of the radial parts of the corresponding single- 
particle wave functions. We then have 

e (N)=e + Nf oa , e+ (N) = e + D + Nf+ a , (35) 

so that the distance D(N, 0) varies as 

D(N,0) = D +N(f+ - f 00 ) ■ (36) 

Significantly, in both the atomic and nuclear problems, 
the sign of the difference /o+ — foo is negative^- so that 
the function D(N, 0) falls off as the number N of added 
neutrons increases. The standard FL scenario, in which 
all added quasiparticles occupy the level 0, remains valid 
as long as the level distance D(N, 0) remains positive. As 
seen from Eq. (|36p. this distance changes sign when N 
reaches the critical number N c — D /(foo — f+o)- Thus, 
at N > N c , 

D(N, 0) = (N - N c ) (f +0 - f 0Q ) < , (37) 

forcing all the quasiparticles to resettle into the level +. 

To perform the next iteration, one calculates the dis- 
tance D(0, N) = e (0, N) — e+(0, N) between levels and 
+ for the case that all the quasiparticles occupy level +, 
obtaining 

D(0, N) = D +N (f+o - f++) ~ -(N-N c ) (f 00 - f+ ) . 

(38) 

Since N exceeds the critical number N c , the sign of this 
quantity is negative. This implies that the priority for 
filling reverses again, requiring that all the quasiparticles 
return to the level 0. This 2-cycle is repeated indefinitely. 

In the weak coupling limit where foo — ► 0, the critical 
number N c diverges, and filling of the level is completed 
before N reaches N c . In this case, level-filling occurs nor- 
mally, the 2-cycle being unattainable. However, in the 
atomic problem, the magnitude of the matrix elements 
of the Coulomb interaction between the orbiting elec- 
trons significantly exceeds the distance between neigh- 
boring single-particle levels. As a consequence, 2-cycles 
emerge at N ~ 1 in the iterative maps for atoms of al- 
most all elements not belonging to the principal groups. 
The implication is that the electronic systems of these 
elements do not obey standard FL theory. The same con- 
clusion is valid for many heavy atomic nuclei with open 
shells. 

Naturally, the occurrence of such 2-cycles in the anal- 
ysis of level-filling in finite Fermi systems presents a 
dilemma that must be resolved. Here (and in the other 
examples) such cyclic behavior is obviously unphysical. 
Here the resolution lies in the merging of single-particle 
levels and their partial occupation, hence in behavior 
that conflicts with standard FL theory yet in fact main- 
tains consistency within the broader framework of Lan- 
dau theory42 
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FIG. 5: Iterative maps for the discretized Nozieres model 
(|39[) when inputs taken from the two preceding iterations 
are mixed with the parameter £ = 0.2. Left panels: spec- 
tra e^'(p) with j — 0,1,2,3, in units of e%. Right panels: 
momentum distributions n^'(p). 



D. Modified Poincare mapping and new insight 
from chaos theory 

Apparently, the occurrence of persistent 2-cycles in the 
iterative maps of Eq. (I2ip prevents us from finding true 
solutions of the fundamental Landau equation ([5]) beyond 
the QCP for a specific class of Fermi systems possessing 
singular effective interactions. It can be argued, however, 
that this failure is a consequence of the inadequacy of the 
iterative procedure employed, which works perfectly on 
the FL side of the QCP. Indeed, a refined procedure that 
mixes iterations does allow one to avoid the emergence 
of these 2-cycles. Nevertheless, the improved procedure 
possesses the same feature: the iterations do not con- 
verge, although the pattern of their evolution becomes 
more complicated and — as will be seen — both intriguing 
and suggestive. 

By way of illustration let us consider a modified 
Poincare mapping for the Nozieres model, again with 
coupling constant g = 0.2e^. Choosing a mixing pa- 
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rameter £, the equation 

e u+1) (p) + = (1 - C) [e (j) (p) + V {J) \ 

+([p 2 /2M + gn {: > ) {p)} (39) 

is used to generate the iterative maps shown in Fig. [5] 
Recovery of the ordinary Fermi distribution uf(p), an 
inherent feature of the standard iteration procedure at 
the second iteration (see Fig. [5]), no longer occurs; in- 
deed, the sequence of iterations fails to converge, even to 
a limit cycle. The number of sheets remains three at the 
second iteration. At the third, however, seven sheets of 
the Fermi surface emerge, and the number of sheets con- 
tinues to increase in successive iterations. At the same 
time, the distance between the sheets continues to nar- 
row, since the domain of momentum space in which the 
improved iteration procedure does not converge remains 
almost the same as that in which the standard procedure 
finds a 2-cycle. This phenomenon of proliferating Fermi 
sheets recapitulates a scenario envisioned in the seminal 
study by Pethick et ali2£ of the quark-gluon plasma based 
on the interaction function (13U1) . 

Treating the number j — 1,2,3,... of the iteration 
as a discrete time step, the sequence of the pictures in 
the left column of Fig. [5] shows the "temporal" evolution 
of the quasiparticle spectrum. At any time-step t, the 
single-particle energy e(p,t) falls off steadily as p goes 
to zero, while its sign changes unpredictably in a finite 
region fl t of momentum space adjacent to the Fermi sur- 
face. These erratic changes of sign induce unpredictable 
jumps of the occupation numbers n(p, t) between the two 
values and 1. Significantly, at T = and as t goes 
to infinity, the region of momentum space in which it- 
erations of Eq. ([5]) fail to converge tends to a definite 
limit fl. The lack of convergence of the iteration process 
may be attributed to the presence of a kind of "quantum 
chaos." Since entropy is a natural byproduct of chaos, we 
tentatively identify 2£lln2 as its measure in the present 
context, where O now denotes the numerical volume, in 
momentum space, of the domain of non-convergence, and 
the factor 2 comes from the two spin projections. Thus, 
2f2 In 2 is interpreted as a special entropy associated with 
a Fermi system for which iteration of Eq. ^ does not 
converge to a solution. 

This naive formula for the special entropy can be re- 
fined. To do so, we introduce a "time-averaged" single- 
particle energy e(p) by appropriating a standard formula 
from statistical physics, 

T N 

£(P) = T l™ \ I e( P , t)dt ee Urn ^ eW) (P) ' ( 4 °) 
{ o 

in which the iterates for e(p) are averaged over the ficti- 
tious time t. In the same way, we introduce a correspond- 
ing time average n(p) of the iterates of the momentum 
distribution n(p). The relation between the two time av- 
erages, simply 

e(p) = p 2 /2M -n + gn(p) , (41) 



stems from Eq. (|33| . Remarkably, the mixing parameter 
C cancels out in arriving at this relation. 

Obviously, n(p) only takes values or 1 wherever the 
iterations converge, and then e(p) is a parabolic func- 
tion of p coinciding with the true single-particle energy. 
However, as seen from Eq. f|40[) and verified by results 
shown in the left column of Fig. [5J the function e(p) van- 
ishes identically in the domain Q where iterations fail to 
converge, yielding 

nM^—^., Pen. (42) 

9 

The notation n»(p) has been introduced to signify a 
smoothed momentum distribution, determined by aver- 
aging iterates e(p, t) for the single-particle energy accord- 
ing to the prescription (|4"0")) . 

Thus, by implementing (i) a modified Poincare map- 
ping procedure based on Eq. ([39]) together with (ii) 
"time-averaging" of iterates for e(p) and n(p) in the man- 
ner of Eq. (|4D|) , we have found a solution of the problem 
that satisfies the Pauli principle in the domain f2 where 
the sequence of iterates does not converge. The solution 
so obtained is independent of the parameters specifying 
the refined iteration procedure. The boundaries pi and 
Pf of the momentum interval pi < pf < Pf defining the 
domain of non-convergence are determined by the condi- 
tions n*{pi) = 1 and ra*(p/) = 0. Thus, the quasiparticle 
momentum distribution n(p) corresponding to the new 
solution, hereafter written as ra*(p), is given by 1 and 0, 
respectively, at p < pi and p > pf, and by Eq. (|42[) in 
between. 

Now we are prepared to introduce a Kolmogorov-like 
entropy 

5* = — 2 J [n*(p) hin*(p) + (1 — n*(p)) ln(l — n lr (p))]dv , 

(43) 

in which the second term of the integrand accounts for 
the equal status of particles and holes. According to the 
definition of n*(p), the integrand in Eq. (|43|) vanishes 
outside the domain Q, and we see that S* is essentially 
proportional to the volume of this domain (also called 
fi). Setting n*(p) inside Q equal to a typical value 1/2, 
one does indeed arrive at = 2fHn2. 

This program is to be implemented similarly when 
dealing with other Fermi systems for which iteration of 
Eq. © does not converge to a solution. The first step is 
to integrate Eq. ([5]) and find an explicit relation connect- 
ing the spectrum e(p) with the momentum distribution 
n(p). The average quantity e(p) is then constructed by 
means of Eq. (|40p . Since this quantity vanishes identi- 
cally in the domain Q where iteration fails to converge, 
we obtain a closed equation 

c(p,n„) = 0, pef!, (44) 

for finding the smoothed, NFL component of the mo- 
mentum distribution n»(p) and calculating the entropy 
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The possibility of a nonzero value of the entropy at 
T = beyond the QCP is an inherent feature of the 
topological scenario being explored within the original 
Landau framework. The existence of such an entropy 
excess for a system whose interaction function / has a 
long-range component means that the associated ground 
state is statistically degenerated These findings offer a 
new perspective not only on the QCP itself, but also on 
the concept of "quantum chaos" — suggesting a novel role 
for the phenomenon of chaos in quantum many-body sys- 
tems. The underlying connections and their implications 
will be developed in depth in a separate paper. 



IV. FERMION CONDENSATION AS A SPECIES 
OF TOPOLOGICAL PHASE TRANSITION 

An equation identical to (|44|) can be derived in another 
based on the premise^ that the ground-state en- 



45 



way, 

ergy E of the system is a functional of the quasiparticle 
momentum distribution n(j>). For canonical Fermi liq- 
uids, the minimum of E[n] is attained at the boundary 
point npip) — 0(p — pp) of the manifold {n} containing 
all candidates for the function n(p) that satisfy the Pauli 
restriction < n < 1 and particle-number conservation. 
In Ref. l45l. the minimum of the model functional 



Z7 f P 2 i \ j , / f "-(PiMPz) > , 
E = J—n( P ) dv + -J 1 -— rdvidv2 



(45) 



was found. The Fermi step np(p) turns out to be the true 
ground-state quasiparticle momentum distribution only 
if correlations are rather weak. Otherwise a topological 
phase transition occurs at some critical coupling constant 
f c - At / > f c , a new distribution is determined from the 
variational condition 45 



SE 
Sn(p) 



p e fi 



(46) 



For the above model functional, this condition takes the 
explicit form 



2M J 



|Pi - P 



dv\ = /j , p g n 



(47) 



This equation is easily solved for n(p) by exploiting a 
transparent analogy with a system of charged particles 
moving in an external elastic field. The result, yielding 
n *(p) = const. < 1 at p < pf and n lr {p) — other- 
wise, is drastically different from the familiar FL solu- 
tion. Importantly, the same result is obtained when the 
iteration/averaging procedure established in Sec. III.D is 
applied to Eq. (j47|) . For the Nozieres model 41 - studied in 
Subsec. III.D, Eq. (|46|) may also be solved analytically, 
and the solution is again in agreement with the result 
(|42| of the procedure introduced there. 

Prior to discussing this coincidence in more de- 
tail, let us focus on Eq. (|4"6"]l itself. This variational 



condition, studied rather extensively during the last 
decade) 41 i 45 . 46 i 47 i 48 i 49 . 50 i 51 i 52 i 53 i 54 i 55 is generic. To illu- 
minate its nature and conceptual status, we invoke a 
mathematical correspondence of the functional E[n{p)\ 
with the energy functional E[p(r)\ of statistical physics. 
If the interactions are weak, the latter functional attains 
its minimum value at a density p determined by the size 
of the vessel that contains it, and which it fills uniformly. 
In such cases, solutions of the variational problem evi- 
dently describe gases. On the other hand, if the interac- 
tions between the particles are sufficiently strong, there 
arise nontrivial solutions of the variational condition 



S_E[p] 
Sp(r) 



(48) 



that describe liquids, whose density is practically inde- 
pendent of boundary conditions. 

The energy functional E[n] of our quantum many-body 
problem must have two analogous types of solutions, with 
an essential difference: solutions n*(p) of the variational 
condition (|46|) must satisfy the condition < n(p) < 
1. This condition cannot be met in weakly correlated 
Fermi systems, but it can be satisfied in systems with 
sufficiently strong correlations. 

Turning now to the coincidence between solutions of 
Eq. (I46p and Eq. (|4~4"1) , we note that the quasiparticle 
energy is by definition just the derivative of the ground 
state E with respect to the quasiparticle momentum dis- 
tribution n(p), or, referred to the chemical potential, just 
e(p) = SE/5n(p) — /i. Consequently, the variational con- 
dition (|46[) may recast as an equation 



e(p, n* 







pel!, 



(49) 



that is in fact identical to Eq. (|44[) . This equivalence 
allows us to use the same notation n* for the NFL mo- 
mentum distribution n(p) in the two formulations of the 
problem. 

The above analysis and discussion make it clear that 
the most salient of the unorthodox features we have un- 
covered in applying the original Landau quasiparticle for- 
malism to the behavior of strongly correlated Fermi sys- 
tems beyond the QCP is that at T = the single-particle 
spectrum e(p) becomes completely flat over a finite do- 
main f2 of momenta adjacent to the Fermi surface. One 
may envision this phenomenon as a swelling of the Fermi 
surface. Another generic feature, concomitant with flat- 
tening of the single-particle spectrum, is partial occupa- 
tion of single-particle states of given spin, i.e., n(p) is no 
longer restricted to the values and 1, but may take any 
value in the interval [0, 1]. 

This extraordinary behavior is associated with a topo- 
logical phase transition fundamentally different from the 
one which features a proliferating number of sheets of the 
Fermi surface, in that the roots of the equation e(p) = 
now form a continuum instead of a countable set. Since 
the Landau quasiparticles experience no damping, the 
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single-quasiparticle Green function becomes 
l-n*(p) u*(p) 



G(p,e) 



e— iS 



pen, 



(50) 



in the domain Q of vanishing e(p,n) and retains its FL 
form for p outside Q. This particular structure of the 
Green function (|50[) may be characterized by a topological 
charge defined as the integra l 28 ' 46 



N = / G(p,£)d ; G(p,£) 



27ri 



(51) 



where the Green function is considered on the imaginary 
axis of energy e — i£ and the integration is performed 
over a contour in (p,£)-space that embraces the Fermi 
surface. For canonical Fermi liquids and systems with a 
multi-connected Fermi surface, the topological charge Af 
is an integer, whereas for the more exotic states charac- 
terized by a completely flat portion of the spectrum e(p), 
its value is half-odd-integral. 2S i 

Aided by the variational condition (|46[) . one can eluci- 
date what happens when a quasiparticle with momentum 
p € f2 is added to the system, assumed again to be homo- 
geneous. In contrast to what happens with a canonical 
Fermi liquid, the addition of just one quasiparticle now 
induces a rearrangement of the whole distribution func- 
tion n*(p). This implies that the kind of system being 
considered cannot be treated as a gas of interacting quasi- 
particles, even though the original Landau quasiparticle 
concept still applies. 

The set of states for which Eq. |gS} [or Eq.(f55))] is satis- 
fied has been called the fermion condensate (FC)^ while 
the topological phase transition in which the Fermi sur- 
face swells from a line to a surface in 2D, or from a surface 
to a volume in 3D, is otherwise known as fermion con- 
densation. In finite systems the phenomenon of fermion 
condensation exhibits itself as merging of neighboring 
single-particle levels. 42 Unfortunately, the terms fermion 
condensation and fermion condensate have unnecessarily 
promoted controversy. Theorists are condition to think 
that in contrast to bosons, fermions cannot condense, be- 
cause fermions cannot occupy the same quantum state. 
However, in everyday life, "condensation" means sim- 
ply a dramatic increase of density. For example, people 
use the word to describe what occurs when the vapor 
in clouds forms into liquid drops that fall as rain. On 
the other hand, statistical physics deals with occupation 
numbers rather than with wave functions, and the ther- 
modynamic properties of a Bose gas at low temperatures 
are properly evaluated by treating Bose condensation as a 
process in which a macroscopic number of bosons occupy 
the same zero-momentum single particle state, all with 
single-particle energy equal to the chemical potential p. 

Ideally, in experimental measurements of neutron scat- 
tering on liquid 4 He, Bose condensation is reflected in a 
sharp peak in the density of states p(e) = p c 5{e), with a 
prefactor p c that is to be interpreted as the density of the 
Bose condensate. In strongly correlated Fermi systems, 



fermions are also capable of condensation in much the 
same sense: a macroscopic number can have the same 
energy p. If so, it follows from Eq. (f5TJ)) that the zero- 
temperature density of states, which is associated with 
the integral of the imaginary part of the retarded Green 
function Gr(p, e) over momentum space, has the same 
kind of peak p(e) — p c 5(e), where p c is the now the FC 
density. We hasten to add that although a macroscopic 
number of fermions have the same energy, they possess 
different momenta; hence the existence of states contain- 
ing a FC does not contradict the Pauli exclusion princi- 
ple. 

In spite of the similarities between boson and fermion 
condensation, an important difference between the struc- 
ture functions describing them must be noted. In the 
Bose case, a macroscopic number of particles reside in 
condensate for all temperatures lower than the critical 
temperature, with all the bosons having zero momentum 
and energy e(p) — 0. On the other hand, in the Fermi 
case, any elevation of the temperature from zero acts to 
lift the degeneracy of the single-particle spectrum in the 
domain Q. Indeed, since a minute temperature elevation 
does not affect the FC momentum distribution, Eq. ((4]) 
can be inverted to yield 41 - 



e(p,T^ 0) =Tln 



1 - n*(p) 
n*(p) 



p e fi 



(52) 



Thus, the dispersion of the single-particle spectrum of 
systems having a FC turns out to be proportional to tem- 
perature, in contrast to the situation for canonical Fermi 
liquids, where it is independent of T. This distinction can 
be employed for unambiguous delineation of the bound- 
aries of the domain f2 occupied by the FC at finite T. 

In nonhomogeneous systems, finite or infinitely ex- 
tended, appropriate single-particle states are no longer 
eigenstates of momentum p, so Eq. (|49|) must be replaced 
by 



e x = o , A e n 



(53) 



where A is set of quantum numbers specifying the single- 
particle state. At finite temperatures, the FC degeneracy 
is lifted according to the same formula (f5"2"]) as for the 
homogeneous case. 



V. NUMERICAL CALCULATIONS 

In this section, we discuss the results of numerical 
calculations of the single-particle spectra e{p,T) and 
quasiparticle momentum distributions n(p, T) beyond 
the QCP. We compare the temperature evolution of these 
quantities in a system for which the interaction function 
f(k) is singular at k — > 0, with their temperature depen- 
dence in systems for which f(k) has no singularities. 

We consider model quasiparticle interaction functions 
/(&) that depend only on the difference k = p pi. In 
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FIG. 6: Single-particle spectrum e(p) (top panels) and deriva- 
tive de(p)/dp in units of v F = p F /M (bottom panels), as 
calculated for the interaction model (|55f) with parameters 
t/2 = 0.16 and fa = 0.14. These quantities are plotted as 
functions of p/p F at four (line-coded) temperatures expressed 
in units of 10 -2 e'p- and at r s = 6.8 (left column), r s — 6.9 
(middle column), and r s = 7.0 (right column). 



this case one can integrate the relation ([5]) over momen- 
tum and arrive at the integral equation 



El 

M 



A* + 2 Jf(P - Pi) n ( e (Pi)) d vi 



(54) 



for the quasiparticle spectrum e(p), the chemical poten- 
tial \x being determined by the normalization condition 
@. To evaluate the T-dependent spectrum e(p, T) enter- 
ing all thermodynamic and transport integrals, we substi- 
tute the formula into the right-hand side of Eq. ([54]) . 
We employ an iteration procedure that mixes iterations 
for numerical solution of Eq. (|54|) . with a mixing parame- 
ter £ having the same meaning as in Eq. (|39[) . It is worth 
noting that although at low temperature the momentum 
distribution n(p, T) is smoothed somewhat, one meets 
the same difficulties with iterative solution of Eq. (j5~4")l 
for singular interaction functions as in the case of T = 
when the mixing parameter £ is much larger than the 
ratio T/e° F . However as soon as £ becomes comparable 
with this ratio, stable convergence of the iteration proce- 
dure is achieved. In the case of non-singular interactions, 
the requirements to be met by £ are less stringent; in 
particular, £ can be taken as large as some tens of T/e F . 
Extremely low temperatures are not accessible to the nu- 
merical analysis since the mixing parameter Q should then 
be taken so small that the CPU time required to attain 
reasonable accuracy is unreasonably large. 

The numerical calculations were performed on a mo- 
mentum grid with the step size 10 _4 pf- The mixing 
parameter ( = 1CP 4 ~ 10~ 2 was taken to achieve sta- 
ble convergence of iteration procedure. The accuracy of 
the numerical solution, measured by the maximum dis- 
crepancy between the l.h.s. of Eq. ([54")) and its r.h.s., was 



fixed at 10 _7 £p'. The number of iterations necessary to 
reach this accuracy was about 10 2 -=- 10 4 , depending on 
£. Convergence of the iteration procedure is a warranty 
that an obtained solution provides a local minimum of 
the free energy functional F[n] = E[n] — S[n)T — fj,N. 

We first describe the situation for regular interactions. 
If the analysis carried out in Sec. Ill is correct, then at 
T = the condition (f2T))) stays in effect, but the system 
possesses a multi-connected Fermi surface. At low T be- 
low the new temperature scale T m , the thermodynamic 
properties of the system still follow standard Landau FL 
theory. However, the enhanced density of states, whose 
value we have determined to be inversely proportional to 
the difference \p — Poc\i produces a great enhancement of 
key thermodynamic characteristics, notably the spin sus- 
ceptibility x(T -> 0) and the ratio C{T)/T at T -> 0. If 
the temperature T reaches values comparable to T m , the 
sharp kink structure seen in the momentum distribution 
n(p, 0) associated with the multi-connected Fermi surface 
becomes smeared, and at T > T m the function n(jp, T) 
becomes continuous and almost independent of T in the 
momentum interval [pi,P3] between the sheets. Accord- 
ing to Eq. ([5^)1 , the dispersion of the single-particle spec- 
trum e(p,T) then becomes proportional to T, as in sys- 
tems with a FC. All these features are confirmed in the 
calculations. 

In Figs. [SHTIJl we show results from numerical calcula- 
tions of the spectrum e(p,T) based on Eq. J5J), for two 
different interaction functions / that are regular in mo- 
mentum space. The choice of the first of these, 



/(*0 



1 



-92 



M(k 2 /4p 2 F -l 



PI 



(55) 



is motivated by the fact that for the 2D electron gas, 
it can adequately describe the results of microscopic 
calculations 29 of the zero-temperature single-particle 
spectra e(p,T = 0) on the FL side of the corresponding 
QCP, i.e. at r s < r w = 7.0. (Again, the dimensionless 
parameter r s is the radius of the volume per particle mea- 
sured in units of the atomic Bohr radius, thus related to 



the Fermi momentum by r a 
interaction, 



/(*) = .93 



n 2 p F 



V2Me 2 /p F .) The second 



(56) 



M k 2 + p 2 p 2 F ' 



which is relevant to a 3D system, is chosen because it 
was employed in the first work^i addressing the deep con- 
nections between the two types of the topological phase 
transitions considered in the present article. The results 
of calculations with these two interactions are compared 
with those obtained from analytic solution, at T = 0, 
of a model of fcrmion condensation in a 3D having the 
singular interaction function^ 



/(A) = 9: 



9 

7T e 



M 



(57) 



Fig. [S] displays the single-particle spectrum e(p) and 
the group velocity de(p)/dp, as calculated on the FL side 
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FIG. 7: Occupation number n(p) (top panels), single-particle 
spectrum e(p) in units of e° F (middle panels), and ratio e(p)/T 
(bottom panels) for the interaction model (|55|l at r a — 7.1 
(left column) and r 3 = 7.2 (right column), both exceeding the 
QCP value — 7.0. Interaction parameters are the same 
as for Fig. [5] All three quantities are plotted as functions of 
p/pF at different temperatures (measured in units e%) below 
the crossover temperature T m = 10 _3 e^. 



of the QCP for the model corresponding to the inter- 
action (|55|) . The group velocity behaves as a parabolic 
function of momentum p. With increasing r s , the bottom 
of the parabola gradually moves downward, and, when r s 
reaches r OQl makes contact with the horizontal axis ex- 
actly at the Fermi momentum pp — p x . Consequently, 
the spectrum e(p,T — 0) ~ (p — Poo) 3 has an inflec- 
tion point at the Fermi surfaced When evaluated with 
n(p) = np(p), the group velocity changes sign from pos- 
itive to negative as r s passes r Ml in agreement with the 
topological scenario for the QCP.— In the vicinity of the 
QCP, the temperature dependence of the group veloc- 
ity obeys the relation vp(T) ~ T 2 / 3 , again in agreement 
with the result obtained in Ref. 140 . 

Figs. [7] and [8J show numerical results for the quasipar- 
ticle momentum distribution n(p) and the spectrum e(p), 
as calculated for the interaction function (I55|) both be- 
low and above the temperature T m . Inevitably, the value 
of T rn is somewhat uncertain, because the alteration of 
the FL behavior in this temperature region is associated 
with crossover, rather than with some second-order phase 
transition. Nevertheless, comparison of Fig. [7J (T < T m ) 
and Fig. [SJ (T > T m ) reveals striking changes in the 
structure of both of the functions e(p) and n(p). Indeed, 
as seen in the top panel of Fig. [Jj a well-defined multi- 
connected Fermi surface, witnessed by a pronounced gap 



FIG. 8: Same as in Fig. [7] but at temperatures above T m . 

in filling, exists only at extremely low T < 10~ 4 £^, while 
at T ~ T m ~ 10~ 3 £^, the gap in the occupation numbers 
closes. 

On the other hand, upon inspection of the top panel of 
Fig. [51 we observe that within a region region pi < p < pt 
all of the curves for n(p, T) for T > T m collapse into a sin- 
gle one — i.e., the momentum distribution n(p) becomes 
T- independent. This behavior persists until rather high 
T. Fig. [5J also indicates that the range 77 = pf —pi practi- 
cally coincides with the difference-of-roots Ap = p% — p\ 
(cf. Subsec. III. A), so that in accordance with Eq. f28| 
we can write 

V ~ %/poo - P ■ (58) 

Moreover, comparison of the bottom panels of Figs. [7| 
and [8j demonstrates that the huge variations of the ratio 
e(p)/T that are so prominent at T < T m completely dis- 
appear in the FC domain at T > T m . Thus, again we 
discover scaling behavior, within the range p, < p < pj 
all of the curves representing the ratio e(p)/T collapse 
into a single one. 

The same conclusions follow from parallel calculations 
carried out for the interaction (|56|) . For this case, Figs. [9] 
and [TOl dedicated respectively to T < T m and T > T m , 
trace the behavior with temperature of the spectrum e(p) 
and momentum distribution n(p). 

The scaling behavior observed in Figs. [TUTOlis also ev- 
ident in Fig. Ill) where corresponding results of calcula- 
tions based on Eq. ([5]) and the interaction model ([57]) 
are displayed. The single noteworthy difference is that 
the linear-in-T dispersion of the single-particle spectrum 
e(p) characteristic of the FC domain is already in effect 
at T — > 0, i.e., T m = in this case. The results of similar 
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FIG. 9: Occupation numbers n(p) (top panel), single-particle 
spectrum e(p) in units of 10~ 3 e% (middle panel), and ratio 
e(p)/T (bottom panel), plotted versus p/pf at three line-type- 
coded temperatures in units of e%, taken below the transition 
temperature T m — 3 x 10 -3 £%. Interaction model (|56|l is 
assumed with parameters gs = 0.45 and f3s = 0.07. 



numerical calculations for the singular interaction func- 
tion (|5TJ)) . to be published separately, support the same 
general conclusions, in particular with respect to scaling 
behavior. From the evidence gathered for different sin- 
gular interaction functions, we infer that scaling features 
that govern the flattening of single-particle beyond the 
QCP and at T > T m are universal. The specifics of the 
interaction function affect only the value of the crossover 
temperature T m . 



VI. RELEVANCE OF THE QCP TOPOLOGICAL 
SCENARIO TO REALITY 

Here we provide a brief assessment of the relevance of 
the theory outlined in this article to the contemporary 
experimental situation. At the outset we should empha- 
size that this theory is applicable only to systems whose 
properties are aptly described by the conventional FL ap- 
proach on the weakly correlated side of the QCP. If this 
condition is not met, the theory is irrelevant. 

At present, experimental information on the QCP is 
available for three types of strongly correlated Fermi 
systems. The first type is a silicon inversion layer 



FIG. 10: Same as in Fig. [9] but at temperatures above T, 



in MOSFET, with the electrons forming a homoge- 
neous 2D liquid. It is known from measurements of 
the magnetic susceptibility and the Shubnikov-de Haas 
oscillation s 56 ' 57 ! 58 ' 59 that the 2D electron gas does obey 
standard Landau theory on the FL side of the QCP. 
Furthermore, recent measurements^ carried out for two 
different versions of the crystal structure, namely (100) 
and (111), and for differing disorder have nevertheless 
demonstrated that the density at which the effective mass 
M* diverges remains the same, ~ 0.8 x 10 n cm~ 2 . 
This finding implies that electron- electron interactions 
rather than disorder are responsible for the occurrence of 
the QCP in the silicon inversion layer. Moreover, the ex- 
perimental value of the Stoner factor, which specifies the 
proximity to the ferromagnetic phase transition, is shown 
not to be enhanced relative to its value in the volume of 
the silicon crystal^ This finding implies that ferromag- 
netism is irrelevant to the rearrangement of the ground 
state that takes place at the QCP. From the theoreti- 
cal side, we have already mentioned that results 29 from 
microscopic calculations bearing on the divergence of ef- 
fective mass of the 2D electron gas are in agreement with 
the topological scenario for the QCP. All of the above 
factors, both empirical and computational, speak for the 
relevance of the proposed topological scenario to the re- 
ality of the 2D electron gas in the vicinity of the QCP. 

Salient experimental information is also available for a 
second class of Fermi systems having a QCP: films of 3 He 
atoms lying on different substrates. With certain qualifi- 
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10 and g a = 70.0. All three quantities are shown as functions 
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The third and most extensive class of systems possess- 
ing a QCP is that of the heavy-fermion metals, in which 
the crystal lattice gives rise to anisotropy of the Fermi 
surface. We have argued in Subsec. II. B. 2 that short- 
range antiferromagnetic fluctuations cannot be relevant 
to the QCP in these systems. On the other hand, the 
existence of topological phase transitions in this class of 
systems is well known ; 30 ' 31 as pointed out in the same 
subsection. A priori there appear to be no serious obsta- 
cles to implementation of the topological QCP scenario 
for the heavy-fermion metals. 



VII. UNCONVENTIONAL 
THERMODYNAMICS BEYOND THE QCP 

Theoretical results based on the ideas and method we 
have developed allow one to reproduce available experi- 
mental data on properties of several concrete strongly- 
correlated Fermi systems in the immediate vicinity of 
the QCP (for example, see Refs. I4p||53l65h . However, 
such validation is to some extent inconclusive, since as 
a rule there exist at least two alternative models that 
are capable of explaining the data at the same level of 
accuracy. Accordingly, in this section we will focus on 
certain unusual features of the thermodynamic proper- 
ties of strongly correlated systems beyond the QCP that 
emerge in the proposed topological scenario. Detailed 
comparisons of the theory with experimental data is re- 
served for a future article. 



A. Entropy excess 



cations, these systems may be modeled as 2D liquid 3 Hc. 
If the dimensionless Landau parameters specifying the 
effective interaction between 3 He atoms are rather small, 
the laboratory systems may found to obey standard FL 
theory. There is experimental evidence that this is indeed 
the case^i&i — in spite of the violation of the Galilean in- 
variance and the presence of disorder, effects introduced 
by the substrate. As reflections of a divergent density 
of states, the presence of the QCP in 3 He films is ex- 
hibited in measurements of the spin susceptibility Xr^ 
and the specific heat with the quantitative results 

showing a substantial dependence on the nature of the 
substrate. A crucial feature of the measurements^ is 
the curious Curie-like behavior^ found for the spin sus- 
ceptibility: x{T,p) = Cos(p)/T, with the effective Curie 
constant C e ff depending strongly on the density p in the 
QCP region. The challenge of explaining this experi- 
mentally documented crossover from Pauli-likc to Curie- 
like behavior is still open. A microscopic understanding 
is necessarily complicated by interactions between sub- 
strate and helium atoms and by layering effects at higher 
areal densities. In the next section we will see that the 
topological scenario for the QCP shows promise of cap- 
turing the essential physics of the crossover. 



In standard FL theory where the entropy S is given 
by Eq. ([3]), the curve S(T) starts at the origin and be- 
gins to rise linearly with T. The situation changes at 
the QCP. In the conventional scenario, which is pred- 
icated on a vanishing quasiparticle z-factor, one has 
S(T) ~ C(T) ~ rin(l/T). In the topological scenario 
for the QCP, the divergence of the ratio C(T)/T is even 
stronger^ . C(T)/T ~ T~ 2 / 3 . The standard FL behav- 
ior of C(T)/T does re-emerge beyond the QCP, where 
the density of states is again finite. However, recovery of 
this behavior occurs only at extremely low temperatures, 
with the proviso that T < T m . At higher T, the FC forms 
in the domain f2, and its characteristic NFL momentum 
distribution, given by ro* (jp) , leads to a drastic change in 
the behavior of the entropy S(T)^^^. The basic en- 
tropy formula of the original quasiparticle formalism 
remains intact, but due to the NFL component in n»(p), 
the system is seen to possess a T -independent entropy 
excess S* (p) . Its value does not depend of the manner 
of its evaluation, since Eq. (|43[) and (J3]) provide the same 
result. 

The situation we now face — with the strongly corre- 
lated fermion system having a finite value 5 1 , of the en- 
tropy at T = — resembles that encountered in a sys- 
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tern of localized spins. In the spin system, the entropy 
referred to one spin is simply In 2, while in the system 
having a FC, we have S*/N ~ 77 In 2, where rj — Pf/p is 
the dimensionless FC parameter. 

Numerical calculations demonstrate that within the 
domain O, the momentum distribution n*(p) changes 
rapidly under variation of the total density p. The cor- 
responding nonzero value of the derivative dS*/dp pro- 
duces a huge enhancement of the thermal expansion co- 
efficient (3 ~ dS(T)/dp ~ 77 with respect to its FL value, 
proportional to Consonant with this result, experi- 
mental data 62 show that at low T in many heavy-fermion 
metals, (3 is indeed almost temperature-independent and 
exceeds typical values for ordinary metals by a factor 
10 3 -10 4 . To our knowledge, no theory has previously 
been advanced to explain this enhancement. 



B. Curie- Weiss behavior of the spin susceptibility 

Another peculiar feature of strongly correlated Fermi 
systems beyond the QCP involves the temperature de- 
pendence of the spin susceptibility x{T) = Xo(T)/(l + 
<7 n (T)), where 



Xo = P 2 MT) = -2p 2 e J 



dn(p, T) 
de(p) 



dv , 



(59) 



and go is the spin-spin component of the interaction func- 
tion (which remains unchanged through the critical den- 
sity region). As mentioned before, beyond the QCP the 
standard Pauli behavior of x{T) prevails only at T < T m , 
and its value, proportional to the zero-temperature den- 
sity of states N(0) ~ \p — p^x, turns out to be greatly 
enhanced. 

At T > T m , insertion of n*(p) into Eq. ((59)) yields the 
Curie-like term 



X*(T) = p 



2 Ceff(p) 



T 



(60) 



with an effective Curie constant 



C cS (p) = 2 J n.(p)(l - n*{p))dv (61) 

that depends dramatically on the density 54 ' 63 Since C e ff 
is proportional to the FC parameter n, we infer that 



CefF — . 



(62) 



Thus, all compounds in which the spin susceptibility ex- 
hibits the Curie-like behavior possess a large entropy. 
Furthermore, in the whole temperature interval from 
T = to T > T m , the spin susceptibility of a Fermi 
system beyond the QCP possesses Curie- Weiss-like be- 
havior x(T) ~ l/(r — TV) with a negative Weiss tem- 
perature TV. Measurements in 3 He films on various 
substrates^ 4 - and in numerous heavy-fermion compounds 
(cf. Refs. Illllolh provide examples of this NFL behavior. 



We emphasize that in our scenario, the negative sign of 
TV holds independently of the character of the spin-spin 
interaction, repulsive or attractive. (In the latter case, 
the magnitude of this interaction must not exceed cer- 
tain limits, as indicated below in Subsec. VIII. B.) At 
the same time, in a system of localized spins the Weiss 
temperature has a negative sign only if the spin-spin in- 
teraction is repulsive. In this case, however, the Stoner 
factor must be suppressed. The necessity of reconciling 
the negative sign of TV with the enhanced Stoner factor 
observed experimentally in the vicinity of the QCP 8,9 
creates insurmountable difficulties for the standard col- 
lective QCP scenario. 

Another conspicuous feature of the physics beyond 
the QCP is associated with the Sommerfeld- Wilson ratio 
R sw = Tx(T)/i4C(T) ~ S*/C{T). Since the excess en- 
tropy S* does not depend on T, it makes no contribution 
to the specific heat C(T); consequently one should see a 
great enhancement of Rsw- 



VIII. RELEASE OF ENTROPY STORED IN 
THE FERMION CONDENSATE 



The diversity of phase transitions occurring at low tem- 
peratures is one of the most spectacular features of the 
physics of many heavy-fermion compounds. Within the 
standard collective scenario ) 16 ' 17 it is hard to understand 
why these transitions are so different from one another 
and their critical temperatures are so extremely small. 
However, such diversity is endemic to systems with a FC. 
Its source may be traced to an obvious fact: The exis- 
tence of the excess entropy 5* at T — would contradict 
the third law of thermodynamics (the Nernst theorem). 
We may recall that in order to relieve themselves of ex- 
cess of the entropy, systems of localized spins order mag- 
netically due to spin-spin interactions. The situation in 
systems with a FC is similar, but there are many ways to 
release the excess entropy 5» as T goes down to zero. One 
possible route for eliminating 5*, already considered, is 
the crossover between the state with a FC and a state 
with a multi-connected Fermi surface as the temperature 
drops below T m . However, this is by no means a unique 
prescription for removing the entropy excess. Other op- 
tions are associated with second-order phase transitions, 
involving violation of a symmetry of the ground state. 



A. Superconducting phase transitions. 

It is instructive to begin with superconducting phase 
transitions, considered already in the first article 4 ^ de- 
voted to fermion condensation. A necessary condition 
for a superconducting transition to come into play is that 
its critical temperature T c exceeds T m . To determine T c , 
one sets — > in the well-known BCS gap equation, 
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yielding 

A L (p,T 0) = / V ^-^^ A^T^ 

(63) 

where Vl is the effective interaction between quasipar- 
ticles in the L-wave pairing channel. It can be demon- 
strated that the FC contribution to the integrand on the 
right-hand side of this equation is dominant. After some 
algebra employing the identity 

tanh(e(p)/2T) _ 1 - 2n*(p) 

e(p) " rin[(l-n,(p))/n,(p)] ' PG ' 

(64) 

one obtains 

, f 1 - 2n*(p) . , 

T c = V L (p F ,p F ) r --- — jr^-rdv . (65) 

21n[(l - n*(p))/»i*(p)] 

Confining the integration to the small FC domain, one 
arrive at— 

T c ~ X V e F , (66) 

where A is a dimensionless pairing constant and 77 ~ 
\p — Poo] 1 ^ 2 is the FC parameter. We see now that 
a remarkable situation arises if the pairing constant is 
large enough to ensure satisfaction of T c > T m . In con- 
trast to the exponentially small BCS critical temperature 
T c ~ e~ 2 / A , the critical temperature in a system having 
a FC turns out to be a linear function of A. 



B. Phase transitions in the particle- hole channel 

Along the same lines, we can consider the possibility 
of a collapse of particle-hole collective modes in systems 
with a FC. Here we restrict ourselves to long- wave tran- 
sitions that give rise to a deformation of the Fermi sur- 
face. Manipulations similar those previously applied lead 
to the relation 

Ol(p) = - J ^P^^^OMdvt , (67) 

where Ol{p) is shape function characterizing the de- 
formation. In the FC region we have dn(p) / de(p) = 
n»(p)(l — n»(p))/T, so that upon retaining only the FC 
contribution Eq. (|67j) is recast as an equation for the 
transition temperature, 

T N = -T L (p F ,p F ) I n*(p)(l - n*(p))dv ~ -fLV^F ■ 
Jn 

(68) 

On observing that the dimensionless Landau parameters 
Fl = fh/N(Q) keep their values through the critical re- 
gion, we may express the transition temperature in the 
form T N = F L r]/N(p). Remembering that the FC den- 
sity r\ is proportional to \p — Pool 1 ^ 2 , while the density 
of states behaves as N(0) ~ \p — Pool -1 , the estimated 



value of Tjv appears to be proportional to \p — Poc\ 3 ^ 2 , 
being comparable with the crossover temperature T m . 
In the case of ferromagnetism, the condition T/v > T m 
is met only if the spin-spin component of the interac- 
tion function is negative and sufficiently large. Other- 
wise, the system avoids the ferromagnetic phase transi- 
tion and the spin susceptibility obeys the Curie- Weiss law 
x(T) ~ 1/(T — TV) with a negative Weiss temperature 
T w . 

Remarkably, in either of the above candidates posed as 
a mechanism for release of the excess entropy 5* through 
a symmetry-breaking second-order phase transition, the 
unconventional state with a FC always corresponds to the 
high-temperature phase, while the low-temperature phase 
possesses more familiar properties. Phase transitions oc- 
cur in any channel where the sign of the effective inter- 
action is suitable for the transition, provided the value 
of the effective coupling constant is sufficient to produce 
the inequality T c > T m . Otherwise, the entropy excess is 
released through the crossover leading to formation of a 
multi-connected Fermi surface. We see, then, that the di- 
versity of phase transitions in the topological scenario for 
the QCP is due to the accumulation by the FC of a big 
entropy at extremely low temperatures. Because of the 
smallness of the parameter \p — poa\, the transition tem- 
peratures appear to be very low and the phase transitions 
themselves inevitably of quantum origin. Several phases 
specified by different order parameters can in principle 
coexist with each other, giving rise to states of an intri- 
cate nature. As the temperature decreases to zero, these 
phases can replace one another. 



IX. CONCLUSION 

A basic postulate of standard Fermi liquid theory 
reads: "At the transition from a Fermi gas to a Fermi 
liquid, the classification of energy levels remains un- 
changed" (Ref. 23, p. 11). Some 20 years ago, the first 
paper on the theory of fermion condensation^ demon- 
strated that this postulate is incorrect beyond a quantum 
critical point, while maintaining the essence of the quasi- 
particle picture. The paper triggered a wave of criticism 
and disbelief; the judgment, "This theory is an artifact of 
the Hartree-Fock method," was typical. By now, debates 
on the subject have become pointless: numerical calcula- 
tions based on Eq. , carried out during the last decade 
and discussed in the present work, provide the best way 
to answer critics. 

It has been a principal goal of this article to investi- 
gate the structure of the Fermi surface beyond the quan- 
tum critical point within the original Landau quasipar- 
ticle pattern^ We have shown that at T = there are 
two different realizations of this pattern, both of which 
involve a topological phase transition. The first is as- 
sociated with the emergence of a multi-connected Fermi 
surface. The second entails the formation of a fermion 
condensate, implicated by the emergence of completely 
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flat portion of the single-particle spectrum e(p) that may 
be envisioned as a virtual swelling of the Fermi surface. 
Such an inflation of the Fermi surface can occur if the 
Landau interaction function / contains components of 
long range in coordinate space. 

We have performed a series of numerical calculations 
that serve to demonstrate the existence of a crossover be- 
tween the two types of topological structure at extremely 
low temperature, an effect that introduces a new energy 
scale T m into the problem. At T > T m , both the momen- 
tum distribution n(ja) and the single-particle spectrum 
e(p) exhibit universal features inherent in strongly corre- 
lated Fermi systems beyond the QCP. To be specific, scal- 
ing features pertaining to the behavior of the momentum 
distribution n(p) and the (scaled) single-particle energy 
e(p)/T in the vicinity of the Fermi surface are indepen- 
dent of the assumed form of the interaction function / 
(and notably whether on not it contains long-range com- 
ponents). The choice of / does, however, affect the value 
of T m . 

We have explored intriguing features of Poincare map- 
ping as a technique for iterative solution of the nonlinear 
integral equation (|SJ) that connects the group velocity and 
the quasiparticle momentum distribution at zero temper- 
ature. For regular interaction functions, the sequence of 
iterative maps converges to the true solution beyond the 
quantum critical point and describe a new ground state 
with a multi-connected Fermi surface. However, if the 
Landau interaction function / has a component of long 
range in coordinate space, 2-cycles are generated in the 
standard iteration procedure, a behavior symptomatic of 
the inadequacy of this solution algorithm as well as the 
failure of standard Fermi liquid theory beyond the quan- 
tum critical point. Adopting a refined iteration proce- 
dure (in which the input for the next step is a mixture of 
the outputs of the two preceding iterations), the spuri- 
ous 2-cycles no longer arise. However, the new sequence 
of iterative maps acquires chaotic features in a region 
of momentum space adjacent to the Fermi surface that 
coincides with the domain in which the 2-cycles were pre- 
viously found. 



We have elaborated a special procedure for averag- 
ing the successive outputs of the iteration process and 
demonstrated that the averaged single-particle energies 
and occupation numbers so obtained coincide with those 
inherent in states with a fermion condensate. The excep- 
tional states involving a fermion condensate are shown to 
possess a nonzero entropy at zero temperature. This re- 
sult does not contradict basic laws of statistical physics if 
and only if the ground state is degenerate. In effect, the 
comprehensive analysis of Eq. ([5]) presented here affirms 
the consistency of the properties of states possessing a 
fermion condensate. The new insight into the role of 
chaos theory in the quantum many-body problem that 
has been gained through our analysis lays a basis for fu- 
ture studies with the potential of wide implications. 

Finally, we have investigated pathways for releasing 
the entropy excess stored in the fermion condensate that 
are associated with different quantum phase transitions — 
necessarily of quantum origin since the transitions in- 
volved occur at extremely low temperatures. 
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